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ABSTRACT 


SEA PLOT is a three-dimensional graphics program written 
for the Compass Hammer research effort. The program portrays 
a War-at-sea scenario in which a single aircraft attacks a 
Single ship target located at the origin of an earth-fixed 
coordinate system. To plot the three-dimensional scenario 
on a two-dimensional plotting surface, a conic projection 
technique is used in which a cone of vision is generated 
based on the position of a viewer's eye. The position of 
the eye 1s automatically placed based upon the initial 
coordinate position of the attacking aircraft. The program 
1s written using VERSATEC software available on the NPS 
IBM 3033 computer. However, the program has been structured 
in such a manner as to permit easy transition to a more 
sophisticated, picture oriented, graphics language such 
as PLOT 80 which would permit an interactive graphics 
employment of the program and possibly animation of the 


attacking scenario. 
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ieee PRODUCT ION 


This thesis provides a program description and user's 
manual for the graphics package SEA PLOT designed for use 
in the Compass Hammer research effort. SEA PLOT is a 
three-dimensional graphics program which provides the user 
a perspective display of a war-at-sea scenario in which 
a single ship target is attacked by a single aircraft. 

The plot depicts (1) ship position and heading, (2) air- 
@eagec elveht profile including coordinate position, roll, 
pitch and yaw, and (3) ship defensive missile launch envelope. 
All information is displayed in a three-dimensional format 
uSing a conic projection in which an optimal viewing 

position is automatically selected based upon the attacking 
aircraft's initial coordinate position. Because the plotting 
package is intended to provide a visual overview of the 
attacking scenario, it is not a reliable source for making 
measurements of ranges, bearings, etc. This results from 

the conic projection technique used which displays solid 
objects in three-space much as they would be viewed by the 
Viewer's eye. Consequently, range marks and rings located 

in close proximity to the viewing position will show a 

Wider dispersion than those at a more distant location. 

The conic projection itself is discussed in detail in 


enapter I. 





Plotted output from SEA PLOT consists of range rings 





: from the center of the target, aircraft flight profile, 
target ship orientation, and missile launch envelope. 

Range rings are displayed at 5000-meter increments from 
mmemeemcer Of the target located at the plot center. 
Maximum range is calculated by examination of the attacking 
aircraft's flight profile such that maximum plot range 
corresponds to the maximum separation between target and 
attacker. The target ship is shown in planform only, but 


1s oriented along a user input heading. Three-dimensional 


display of the target ship, which is possible with only 





| 
| 
minor program changes, resulted in a cluttered plot output 
| and consequently was omitted. Display of a defensive 
missile launch envelope consists of a cylinder comprised 
of unconnected dots corresponding to maximum missile slant 
range. The cylinder is truncated at a height of twice the 
maximum aircraft flight profile altitude. Aircraft ground 
track is depicted via a series of dark dots plotted on the 
surface corresponding to aircraft position over ground. 
The air track is shown in perspective view as a bold line 
tracking the vehicle movement through the air. Along the 
air track, an aircraft body is positioned and drawn 
corresponding to a minimum of 2000 meters of ground travel. 
This body plot displays aircraft maneuvers in roll, pitch, 


and yaw. Aircraft altitude perspective is enhanced via 


vertical dashed lines connecting ground and air tracks. 





sen PLOT also contains a variety of input scaling 
parameters which permit the user to focus on the close-in 
encounter between the target and aircraft, or to expand 
the view and observe the entire profile from initial to 
final aircraft positions. This is accomplished by entry 
of the scale parameter, one of the program data input 
elements. For example, SCALE=1 causes a wide angle display 
of the scenario. SCALE=2 in essence magnifies the plot size 
by a factor of two, focusing on the area from the origin 
out to one half the original range. Figures 1-1, 1-2, and 
1-3 contain examples of plot output using SCALE=1, 2 and 
Be respectively, Chapter III discusses nee and the 
required program inputs for SEA PLOT. 

Because SEA PLOT output is to be included in the Compass 
Hammer report, hard copy, printed output was required, 
as opposed to an interactive graphic output onto a CRT 
terminal. Consequently, the plotting package was written 
using VERSATEC software available on the NPS computer. 
Because this software system is oriented primarily toward 
display of data in graphic or tabular form, and not toward 
the drawing of pictures, the programming algorithms were 
more complicated than would be required for an interactive 
graphics software package, such as PLOT 80 used on many 
Hektronix Machines. However, because the program does have 
many interactive graphics applications, such as an 


iterative design of an aircraft mission profile, it has 





been structured as much as possible to permit a relatively 
simple conversion to other software systems. To accomplish 
mis, the graphics instructions are written in internal 
Subroutines, each of which contains the commands peculiar to 
the particular software system employed. For example, to 
draw a line connecting an array of coordinate points, 
SEA PLOT calls an internal subroutine DRAW. Resident in 
DRAW are the VERSATEC commands necessary to draw the line. 
To convert to another graphics system, one would merely 
have to adjust the instructions in DRAW to facilitate the 
software conversion rather than search through the entire 
program locating specialized line drawing instructions. 
Since SEA PLOT was written in as general a form as 
possible, its applications are not limited to the Compass 
Hammer project. With generality in mind, it was prepared 
as a stand-alone program and was not incorporated into one 
of the Compass Hammer computing routines. SEA PLOT, in 
essence, can be utilized by any user desiring a three- 
dimensional attack scenario plot. It need not be constrained 
to the display of attacking aircraft encounters, but could be 
used from the standpoint of a ship defending itself from 
a Single aircraft, evaluation of defensive missile design, 
or evaluation of weapon scenario development. Another 
possible use of SEA PLOT would be in conjunction with the 
POO] Input Program (PIP), which currently only provides 


two-dimensional plotted output. 
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Detailed discussion of SEA PLOT techniques and use 
| follow in Chapters II and II which cover the conic pro- 
jection and program input and organization, respectively. 
Appendix A provides a complete program listing, Appendix B 
contains a listing of the conic projection routines, and 
Appendix C discusses use of the Tektronix developed inter- 
active graphics software PLOT-80 on the Tektronix 4081 


meaphiecs terminal. 
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ieee eer PROJECTION 


This chapter discusses the projection technique employed 
by SEA PLOT to provide a perspective two-dimensional 
display of an object in three-dimensional space. The conic 
projection technique was selected for SEA PLOT from among 
a multitude of algorithms [Ref. 1]. This technique is 
particularly well suited to the purpose of displaying a 
three-dimensional object as if it were being physically 
viewed by the user's eye in that not only is the solid 
body orientation in space accurately portrayed, but per- 
spective depth and distance from the viewer's eye are 
preserved. Thus, it becomes an excellent tool for viewing 
an aircraft-ship engagement as observed from a distant 
hecacion. 

As with any effort to display a solid object using a 
mathematical definition of body doorcinates, the elimination 
of hidden lines and surfaces as they are cloaked from view 
by surfaces in a more prominent viewing position presents a 
problem to the programmer. Moreover, the computation time 
required to mask hidden surfaces grows with the square of 
the situation complexity [Ref. 2]. However, the use of 
Simplistic solid body models as opposed to detailed body 
structures eliminates the need for a comprehensive hidden 


surface removal routine. Since SEA PLOT is intended to 


i 





provide a symbolic overview rather than a definitive 
display of an attack scenario, this technique has been 
adopted. Consequently, the aircraft body is displayed in 
triangular form approximating the wing planform of a delta 
Wing aircraft. The triangle is elongated with an included 
tail section to enhance the perspective orientation of the 
aircraft in three-space. Likewise, the target ship is 
projected as a two-dimensional hull form oriented along 
auser's input target heading. However, should future use 
of SEA PLOT dictate more comprehensive target and attacker 
models, Reference 3 characterizes ten hidden surface algo- 
rithms which a programmer may employ. 

iiesconic projection technique involves placement of 
the viewer's "teye’' in an earth-fixed coordinate system. 
Associated with the positioning of the eye is a cone of 
vision encompassing the field of view and a projection 
plane onto which the scenario is to be displayed. SEA 
PLOT automatically determines these positions based upon 
the input flight profile of the attacking aircraft. The 
viewing plane is placed 2000 meters beyond (as viewed from 
the origin) the attacker's initial flight path point, 
oriented perpendicular to a line drawn through the earth- 
fixed origin (the location of the ship) and the initial 
attacker position. The eye coordinates are placed 2000 
meters beyond the viewing plane along the same line. This 


automatic placement of the viewing system eliminates the 
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requirement for the user to determine an optimal viewing 
point which, if ill-placed, could result in plot distortion 
should the field of view exceed the cone of vision. Should 
another viewing position relative to the attacker's air 
track be desired, Chapter III discusses the procedure to 
alter the perspective position. Appendix B contains the 
@eenal routines employed by SEA PLOT to accomplish the conic 
projection and can be used as a "skeleton program" to be 
applied to three-dimensional graphics endeavors of some 
other nature. 

The conical projection method described in this 
chapter is essentially the same as that used in some of the 
displays with the Differential Maneuvering simulator at the 
NASA Langley Research Center [Ref. 4]. It is summarized in 
Reference 5, from which the following derivations were taken. 
Reference 5 also discusses additional applications of the 
technique. 

First consider an aircraft at some distance from the 
center of an earth-fixed coordinate system. The coordinate 
system, with axes labeled Xe, Ye, Ze is right handed, the 
Xe-Ye plane tangential to the earth at the origin, with the 
Z axis pointing up as shown in Figure 2-1. 

The solid body (in this case the aircraft) center of 
gravity (origin of the aircraft body-fixed coordinates) is 
then positioned at Xe, Ye, Ze, with the body orientation in 


three-space given by the Euler angles w, 8, $ corresponding 


Ly 





Bomiedatneg., pitch, amd roll, respectively. The aircraft 
meeeeaxes (Xb, Yb, Zb) are chosen such that Xb points from 
the origin toward the nose, with the Yb axis pointing in 

the direction of the left wing which gives a conventional 
right-handed system with the Zb axis pointing up from the 
aircraft center. Coordinates of the aircraft body are 

meet Xp(i), Yp(i), Zp(i) which are relative to the air- 
craft center of gravity. These body coordinates can then be 
Gee@ietormed to earth-fixed coordinates Xp'(1), Yp'(i), Zp'(1) 


using 


Xehu(d) Xp (i) Xe 
veep (1) ee) lei (4) a |e 
%'pi(1) Zp (1) Ze. 


where [p| is the direction cosine matrix from earth-fixed 


axes to body axes: 





cos 8 cos wp cos 6 cos jp Sin 98 
- cos wt! sin 9 sin 6 cos w~ cos cos § sin wv 
[D] =}+ sin p cos 6 - Sin b! sin @ Sin 6 
- cos wy Sin 8 cos 6 - Sin ! sin @ cos ¢€ cos @ cos ¢ 
+ sin sin 6 - cos p sin 6 


[ 





18 








After transforming the aircraft body coordinates to 
earth-fixed coordinates, the viewing plane must be established 
@meo which the conic projection will be made. In so doing, 
it 1S convenient to establish a new coordinate system 
Mew 1ev, Zev, with the origin at the eye point (Xeye, 

Yeye, Zeye). The X axis of this coordinate system extends 
from the eye point to some view point located on the viewing 
plane. This, then, represents the axis of the projection 
cone. The viewing coordinate system is obtained by two 
rotations of the earth-fixed system. First, a rotation is 
made about the Ze axis by the angle @ev, the angle between 
the Xe axis and the projection of the line connecting the 
eyepoint and the viewpoint onto the Xe-Ye plane. The second 
1S a pitch rotation about the new Y axis by an angle 8ev. 
This orients the cone of vision in space. From Figure 2-1, 


it can be seen that 


1 ,Zview - a 


Sev = tan” ( RID 


where 


R2D = [(Xview - Xeye) * + (Yview - Yeye) “71/2 


Any point in the earth-fixed coordinate system can 


now be expressed in the viewing coordinate system by 


19 





Xev VoD) Xeye 
ew = (Aly p (2) - Yeye 


Zev bic (ly) Zeye 
[ be 4. 


Peewee {A| 1S the previously defined [{D] matrix with 


89 = 6evV 
w wev 
w) oev 


The final step, after transforming the body coordinates 
of the solid object to viewing coordinates is to make a 
two-dimensional projection of the three-dimensional object 
Gmeema tlat picture plane. (For SEA PLOT, this picture 
plane is the plotting paper; but more generally, it would 
be the screen of the graphics display console.) To do this, 
the distance from the viewer (user's eye) to the picture 
plane is defined as the distance XPLANE. (See Figure 2-2.) 
This permits determination of the two dimensional plotting 
coordinates [YSCOPE, ZSCOPE) for the projection of a three- 
dimensional object onto a two-dimensional plotting plane. 
Feemerigure 2-2, it can be seen that the two triangles 


OSP and OP"P' are similar. Therefore, 


XPLANE _ YSCOPE 
Xev Yev 
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OT 


YSCOPE = *&% x XPLANE 


Xev 





Also similar triangles are OPP' and OPP'. Therefore, 
ZSCOPE _ OP 
ZeV OP! 


From the previous set of similar triangles, it follows that 


OP' _ XPLANE 
OPt Xxev 


Hence, 


ZeV 


Zo COUEm Ror * XPLANE 





Appendix B provides the FORTRAN coding to implement 


these derivations. 


st 








Ficure 2-1: 
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Figure 2-2: Transformation from the EV system into the SCOPE system. 
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IfI. PROGRAM INPUTS AND ORGANIZATION 


This chapter discusses SEA PLOT required input data and 
the organization of the program itself. It provides informa- 
tion on the input data and input format as well as detailed 
discussion of the thirteen different routines that comprise 


Sea rt LOT , 


A. PROGRAM INPUTS AND FORMAT 
Program inputs for SEA PLOT are intended to provide 
the user with a wide variety of options in designing the 
desired display output. SEA PLOT inputs fall into two 
categories: plotting parameters and flight path coordinate 
points. A maximum of 200 flight path coordinate points may 
be plotted. Following is a description of input elements 
and format. 
1. Card 1 - Plotting Parameter Input Card 
a. Col 1-3: Enter integer number of flight 
path milestones, right justified. 
b. Col 4: Enter integer value of the type of 
Ship defensive missile envelope. The two types of envelope 
display available are a solid cylinder (Figure 3-1) ora 
bold black range ring drawn on the sea surface (Figure 3-2). 
Integer inputs are: 
1 - range ring. 


0 - solid cylinder. 


24 


Coote so. sontersinteger value indicating if a 
plot of the coordinate axis is desired as follows: 
1 - axis plot is desired (Figure 3-3). 
0 - no axis plot desired. 
d. Col 6-15: Enter decimal value corresponding 


to the desired plot scale. Figures 1-1, 1-2, 1-3 show 


> 


representative scaled plots corresponding to SCALE=1. ,2.,3. 


> 3 


respectively. This input parameter essentially magnifies the 
mlorewith the center of magnification at the origin of 
the earth-fixed coordinate system. No distortion of the 
perspective view results from increased scale factors. 
e. Col 16-25: Enter floating point ship defensive 
missile range in meters. 
f. Col 26-34: Enter floating point value of 
target heading in degrees, counter clockwise from the x axis. 
2. Cards 2 - 201 (One Card per Milestone) 
a. Col 1-10: Floating point value of aircraft 
x coordinate in meters. 
Bb. Col 11-20: Floating point value of aircraft 
y coordinate in meters. 
—eCimeZ ou Lii@ating point value of alrcraft 
z coordinate in meters. 
Peco 10-) Lleating point value of aircraft 
moll in degrees. 
e. Col 41-50: Floating point value of aircraft 


pitch in degrees. 
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Peo sm o0s vOatine point value of aircraft 


heading in degrees. 


B. PROGRAM ORGANIZATION 

SEA PLOT consists of thirteen different routines 
which can be categorized as plot management routines, 
graphics routines, conic projection routines, and plotting 
routines. This organization is intended to provide a 
smooth flow of data from one program element to another as 
well as to present a reasonably coherent program structure 
to the user. Sufficient commentary data has been included 
in the program to facilitate future changes or alterations 
in aS uncomplicated a manner as possible. Program comments, 
coupled with the detailed discussion in this chapter, should 
provide sufficient program background such that program 
modifications to meet a variety of needs could be relatively 
expeditiously made. However, a word of caution 1S necessary 
when altering data or variable elements. Much data is 
passed from routine to routine in the form of Common Blocks. 
While this technique provides an efficient data flow, it 
does detract from the comprehensibility of the coding. 
Consequently, altering data in one subroutine could have 
unpredictable consequences if it impacts on calculations 
taking place in other subroutines. The majority of 
the information passed in common involves input data arrays 


and scaling parameters, some from external sources and some 
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internally generated. Before altering these variables, 
one should thoroughly understand their uses in the various 
algorithms in which they are employed. 
1. Plot Management Routines 

These routines consist of those algorithms required 
to Input data, determine plotting scales, label plots, and 
set up the eye-view coordinate system. Three routines com- 
prise this section; they are MAIN, SEAPLT, and AXIS. 
Each of these is described below, 

a. MAIN Program 

MAIN is straightforward in that its primary 
purpose is to read the input data in the format previously 
described. It also establishes the common blocks "AIR," 
eet" and "INPUT." Variables and their descriptions 
are as follows: 

(1) XAIR, YAIR, ZAIR - Floating point arrays 
of 200 elements each which contain the attacking aircraft 
%, y, z coordinates, respectively. 

(2) ACROLL, ACPTCH, ACHEAD - Floating point 
arrays of 200 elements each which contain the aircraft roll, 
pitch, and heading corresponding to each coordinate position. 

(3) MDISP - Integer input value determining 
the type of missile envelope display to be plotted. 
MDISP=1 causes a single bold range ring to be plotted 


on the surface (z=0). MDISP=0 causes a three-dimensional 


Ad 





=aeweer tO be plotted. [he height of this cylinder 
corresponds to twice the maximum aircraft input altitude. 
The radius of the envelope is the missile slant range 
Which is an input variable. (See Figures 3-1 and 3-2.) 

(4) IAXIS - An integer input value used to 
determine if a plot of the coordinate axes is desired. 
IAXIS=0 indicates no axes are to be drawn. IAXIS=1 
causes the coordinate axes x, y, z to be displayed. 

(5) SCALE - An input floating point value 
used to determine the magnification factor to be applied 
to the plotted output. This permits the user to "zoom in" 
on the center of the plot, magnifying it by the value of 
SCALE. Figures 1-1, 1-2 and 1-3 demonstrate the use of 
this parameter. 

(6) RNGMIS - An input floating point variable 
used to define the maximum range of the ship's defensive 
ideesa ies. 

(7) THEAD - An input floating point variable 
corresponding to the heading of the ship target measured 
in degrees from the axis. 

Printed output from MAIN is a list of input 
variable values and aircraft flight path data. 
b. SEAPLT 

SEAPLT is the driving routine for the entire 

program. It searches for maximum flight path values, 


positions the eye and viewplane coordinates, and sequences 


Zo 





through the plotting routines to cause the graphic output. 
The viewing plane is positioned 2000 meters behind the 
initial flight path position along a three-dimensional line 
from the origin through the “ing PACS position. The center 
of the eye-view coordinate system is placed on this same 
line 2000 meters beyond the viewing plane. Because of the 
low aircraft altitudes employed in current war-at-sea tactics, 
this routine performs an altitude enhancement which adjusts 
aircraft input altitudes such that the maximum plotted 
altitude corresponds to the viewing height. This preserves 
a perspective in which the viewer is always looking along 
or down upon the flight path. In that all of the flight 
path altitudes are increased by this same value, no plot 
distortion occurs, overall perspective is maintained, and 
altitude maneuvers are enhanced. 

Variables used in SEAPLOT are as follows: 

(1) XVIEW,YVIEW,ZVIEW - x, y, z coordinates of 
the center of the viewing plane. 

i) XEME YEN ZEYE - x) y, 2 coordinates (earth- 
fixed) of the eye position. This becomes the origin of 
the eye-view coordinate system. 

(3) XPLANE - Physical distance from the viewer 
to the plotting surface. This variable is particularly 
germaine to plotting on a graphics terminal. Here it is 


fixed at 30 cm. 


Zo 





(2) Dx — Separation of range rings from the plot 
origin. DX is currently fixed at 5000 meters. 

(5) BESTZ - Altitude enhancement factor applied 
to the input aircraft flight profile altitude. 

(6) XMAX - Maximum aircraft x coordinate. 

(7) YMAX - Maximum aircraft y coordinate. 

(8) ZMAX - Maximum aircraft z coordinate. 

(9) THETA - Angle in radians of the line drawn 
through the origin and the first aircraft coordinate position. 

Gey AXIS 

In addition to drawing coordinate axes, subroutine 
AXIS also performs plot labeling and display scaling. Con- 
sequently, this subroutine must be called whether or not a 
plot of the coordinate axes is desired. The only parameter 
required by AXIS is NFLAG. This is the [AXIS argument pre- 
viously described. 

Scaling performed by AXIS is accomplished in such a 
manner as to center the plot on a Versatec 20-inch by 20-inch 
plotting sheet. To do this, AXIS first projects a circle of 
radius corresponding to the length of a vector drawn from 
mpemearth-tixed Origin out to the point given by the coordi- 
hates (XMAX,YMAX). This circle consists of 180 data points. 
After calling the conic projection routines to make the 
projection circle, AXIS searches for the longest vector from 
mae Orlfpin oUt to the periphery of the circle, Since the 


projected circle will appear as an ellipse whose elongation 
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will depend on the viewing angle, this longest vector will 
Vary from view to view. Once this longest projection vector 
is known, AXIS reorigins the plot to the center of the 
meotting sheet at plotter coordinates (10,10), then scales 
@eemtength of the longest projection vector, named SFACT, 
such that its physical size is 10 inches on the plotter. 
This permits the development of an over-all plot scaling 


factor XSCALE where: 


XSCALE = 10.0/SFACT 


Plot magnification can be accomplished to meet the 
user's specifications by multiplying XSCALE by the input 
variable SCALE. This, of course, will result in the ultimate 
computation of many vectors whose length will exceed the 
dimensions of the plotting paper. However, these vectors 
will be clipped at the edge of the plotting sheet and will 
meemcdcetade the plotted output; though a diagnostic, which 
may be ignored, will be generated on the ouptut listing 
containing the flight path data. 

The final function of AXIS is to label the plot. 
The labeling instructions in the program are software 
dependent. Those in AXIS are written specifically for 
VERSATEC software. Hence, careful attention must be paid 
to this routine when adapting it to another graphics system. 

Axis options and labeling are shown in Figure 3-5 


Which shows the axis output. 








2. Graphaics Routines 
The graphics routimes described in this section are 
those elements which accomplish the actual pen moves required 
to draw lines or dots on the plotting sheet. Thus, these 
routines are highly software dependent and represent the 
major changes required to adapt SEA PLOT to another graphics 
Seorem. SEA PLOT requires only three types of graphics 
operations: drawing solid lines, drawing vertical dashed 
lines, and plotting dots. The two subroutine descriptions 
which follow discuss these operations. 
a. Subroutine DRAW 

This subroutine causes an array of projected 
COGrdinate data points to be connected either by a solid 
line or to be plotted as dots on the plotting sheet. Up 
to 200 data points may be plotted. Projected coordinates 
are first scaled to a 20-inch by 20-inch plot size using the 
XSCALE factor developed in AXIS. They are then shifted 
from the earth-fixed coordinate center to the plotting sheet 
coordinate center by adding a value XORIG, or YORIG which 
corresponds to the displacement of the projected earth 
Sepein £rom the plotter origin. DRAW next causes a pen up 
move to be made to the initial array position and then either 
draws a line connecting the points or plots a dot correspond- 


ing to each coordinate position. 





Variables used in DRAW are: 

(1) X,Y - 200 element input arrays of projected 
data to be plotted. 

(2) XPLOT, YPLOT - 200 element arrays containing 
the scaled projection data to be plotted. 

(3) N - Integer value indicating the number of 
coordinate pairs to be plotted. 

(4) ITYPE - Integer value which determines the 
type of plot to be drawn. ITYPE=0 causes dots to be plotted. 
ITYPE=1 causes solid lines to be drawn. 

b. Subroutine VDASH 

Subroutine VDASH causes a vertical dashed line, 
parallel to the z axis, to be drawn between two points whose 
coordinate values are passed as arguments. It should be 
emphasized that only a vertical line is drawn. Should the 
two points to be connected by the line not be vertically 
positioned, the routine will draw a dashed line from the 
first point up to the Y value of the second point. As with 
subroutine DRAW, all data values are scaled to the XSCALE 
factor generated in AXIS. 

Variables used in VDASH are: 


(1) X1,Y1 - Projection coordinates of the lower 


data point. 

(2) X2,Y2 - Projection coordinates of the upper 
data point. 

(3) DY - vertical separation between the data 
points. 








Se) Conic Projection Routines 

The three subroutines that follow are those routines 
used to implement the conic projection technique derived in 
Chapter II. They are also reproduced in Appendix B to 
serve as a ready reference "skeleton program" to be used in 
other graphics efforts. 

a. Subroutine MULT 

MULT performs a matrix multiplication of a 

square 3 x 3 matrix with a column vector, returning the 
resultant column vector to the calling program. The opera- 


tion performed is 


Xnew xX 

Ynew ee ose 4 

Znew Z 
al 


where the vector[X,Y,Z] and the matrix [D] are passed to 
the subroutine as arguments. 
bee oubroutine DIRGOS 

DIRCOS is a rather simplistic routine which 
merely sets up the direction cosine matrices discussed in 
Chapter II. Inputs to the routine are the Euler angies 
psi, theta, and phi corresponding to solid body roll, pitch, 
and yaw in three-space. DIRCOS output is the three-dimensional 
matrix |[D| oT [al employed by the projection. Also available 


as an output is the transpose of [DI. In that the coding used 
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in DIRCOS is a straightforward assignment of variables to 
matrix positions, no further discussion will be attempted 
here . 

eo Souemomeine CONERO 

CONPRO is the primary routine for the conic 
projection scheme implementing the derivations in Chapter II. 
Inputs to CONPRO are the center of gravity coordinates of 
Miemaoll1d body to be projected, the number of defined body 
coordinates, the array of body coordinate points, and the 
body Euler angles in degrees measured counterclockwise from 
the earth-fixed X axis. After performing all translations 
and rotations required by Chapter II, CONPRO outputs the 
plotting surface coordinates of the projected solid body. 
These coordinates are passed by the calling program to the 
DRAW routine where they are scaled and plotted. 

Variables used in CONPRO are: 

(1) CG - Three element array containing the 
X,Y,Z earth-fixed coordinates of the solid body center of 
aravit Yu. 

(2) NPOINT - Integer value indicating the 
number of data points used to define the solid body. 

(3) XP - Array of body X coordinate values in the 
body-fixed reference system. This array is variably dimen- 
Sioned to NPOINT elements. 

(vee ina, Or body Y coordinate values im the 
body-fixed reference system. This array contains NPOINT 


elements. 
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(5) ZP - Array of NPOINT body Z coordinate values 
in the body-fixed reference system. 

(6) EULER - Three element array containing the 
Euler angles for roll, pitch, and heading in degrees. 

(7) YSCOPE - Output array of NPOINT projected 
plotter X coordinates. 

(8) ZSCOPE - Output array of NPOINT projected 
plotter Y coordinates. 

4. Plotting Routines 
The five subroutines that follow are those algorithms 

which develop the different picture segments that comprise 
the SEA PLOT output. The five elements displayed are: the 
circular surface grid consisting of 5000-meter range rings, 
Meemaircrait ground and air tracks, the aircraft body plot, 
the target body plot, and the defensive missile launch 
envelope. Since more sophisticated graphics languages such 
as PLOT 80 employ the concept of graphics control blocks 
(discussed in Appendix C) or picture segments, SEA PLOT 
was organized along these same lines to facilitate easier 
transition to one of these more powerful software systems. 
Using this scheme, each subroutine will cause a specific 
picture segment to be displayed as directed by the sequencing 
routine SEAPLT. As a result, the user has the option of 
deleting those picture elements which are not desired. For 


example, by circumventing the call to the target subroutine, 


36 








one could get a plot of only the surface grid and aircraft 
flight profile without a target being shown, 
a. Subroutine CGRID 

Subroutine CGRID projects and draws a series of 
concentric range circle centered at the target position at the 
Origin. As presently written, these range circles are drawn 
at 5000-meter increments as indicated by the variable DX 
discussed in SEAPLT. Each range circle is defined by plane 
polar coordinates consisting of 181 data points per circle. 
The number of circles drawn is a function of the maximum 
X coordinate for the aircraft flight profile such that the 
number of circles is the quotient of the maximum X range 
and the radius increment DX, 

Variables used in CGRID are: 

(1) X,Y¥,Z - 181 element arrays containing the 
unprojected coordinates of the data points defining each 
range ring. The Z array is set to zero corresponding to 
Zomeomaltitude . 

Pest Ore coeere - Isl-element arrays containing 
the projected plotter coordinates of the range circle data 
points. These arrays are passed to the DRAW subroutine where 
enecy are scaled and plotted. 

(3) DX = Tmterval between range circles in the 
earth-fixed reference system. 

(4) RADIUS - Radius of each range circle in the 


earth-fixed system. 
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Figure 3-4 portrays the output from CGRID. 
DeESUDLoOur Ine AcPLOT 

Suibpeoucine AGPEOP plots up to 200 aircraft 
flight path coordinate points in a three-dimensional 
perspective view. To enhance the view, both ground and air 
tracks are displayed. The ground track is plotted by setting 
micmamrcratt Z coordinate to zero, then projecting the X,Y,Z 
coordinates and plotting the data as series of dots corre- 
Seomarme to each input point. In a like manner, the air 
track iS projected without zeroing the Z coordinate. A 
Semra. Dold line connects all air track data points. To 
memmeorce the altitude perspective, an optional dashed 
vertical line is drawn between corresponding air and ground 
track data points. This vertical line option parameter is 
a program change option which is specified by the calling 
program not by a user input. 

Aircraft maneuvers in roll, pitch, and yaw are 
displayed by plotting the aircraft body at positions corre- 
Sponding to every 2000 meters of ground travel in the earth- 
fixed system. This is accomplished by a check of ground 
distance traveled between input points followed by a call 
to subroutine ACRAFT after 2000 meters of displacement has 
been detected. Because the surface area of the aircraft 
Bedyetseinsigniticant with respect to the total surface area 
aesmlayed by the plot, a magnification factor is applied 


to the aircraft body. It has been experimentally determined 
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EEE, 


that a magnification factor of 15 permits an optimal body 
display. 

Inputs to ACPLOT include the number of milestones 
memee plotted, the vertical line option parameter, and 
Maenitication factor. These and the other internal variables 
used in this routine are discussed as follows: 

(1) N - Integer number of aircraft milestones 
to be plotted. 

(2) VLINE - Decimal flag used to exercise the 
Geemea vertical line option. if VLINE=1, the vertical 
altitude lines are drawn. Otherwise, no vertical lines 
are drawn between tracks. 

(3) AMAG - Decimal magnification factor to be 
applied to the aircraft body plot. 

(4) X,Y,Z - 200 element arrays containing the 
earth-fixed center of gravity coordinates (input milestones) 
for the body plot and the flight profile tracks. 

(5) XPLOTA,YPLOTA - 200 element arrays con- 
taining the projected X,Y coordinates for the air track. 

(6) XPLOTG,YPLOTG - 200 element arrays con- 
aaimine the projeeted X,Y plotter coomdiymates for the ground 
track. 

Figures 3-5 and 3-6 show the plot output options 


neue AGPLOT . 
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Se ne =e 


@  Sibllpacteibi cakaye, ACRAFT 

ACRAFT contains the predefined aircraft body 
coordinates in the body-fixed coordinate frame, This routine 
applies the magnification factor to the body coordinates, 
causes a projection of the body to be made and a plot output 
to be generated. The actual body coordinates are defined 
in data statements in the program and are not user verifiable. 

Inputs to ACRAFT are the Euler of angles roll, 
pitch and yaw, the center of gravity coordinates in the earth- 
fixed reference, and the magnification factor. These and 
other variables used by ACRAFT are defined as follows: 

(1) EULER - A three element array containing 
the body roll, pitch, and yaw angles in degrees. 

(2) CG - A three element array containing the 
X,Y,Z earth-fixed coordinates of the center of gravity. 

Perne, YACW 2AC = Eight ellement arrays 
containing the predefined aircraft body coordinates in the 
body reference frame. 

(4) AMAG - Magnification factor to be applied 
to the body coordinates. 

The output from ACRAFT can be seen in Figures 
3-5 and 3-6. 

d. Subroutine TARGET 

TARGET performs essentially the same function as 

ACRAFT except that it contains the predefined body coordinates 


of the ship target. Because an excessive number of plotted 
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vectors tends to make the overall plot output too busy, the 
target is displayed as a hull planform only. The planform 
meesiaded tO accentuate its position on the plot. At present, 
the target is plotted at the earth-fixed origin. However, 
Gems is not a rigid requirement since this particular 
routine has been written in a general manner to permit 
positioning the target anywhere in the viewing area. Thus, 
the location of the target center of gravity is an input 
from SEAPLT, the calling program, and can be altered to fit 
the user's needs. The other input to TARGET is the ship 
heading measured in degrees from the earth-fixed X axis. 
This also is a user specified input parameter. 

Variables used in TARGET are: 

(1) CG - A three element array containing the 
target earth-fixed center of gravity coordinates. 

(2) THEAD - Target heading measured in degrees 
counterclockwise from the earth-fixed X axis. 

(3) TMAG - Magnification factor applied to the 
body fixed target coordinates. This variable is set to the 
same value as AMAG in ACRAFT. 

(4) XTGT,YTGT,ZTGT - 26 element arrays containing 
the predefined body coordinates of the target in the body- 
fixed reference system. 

(5) XSCOPE,YSCOPE - 26 element arrays containing 
the projected coordinates of the target in the plotter refer- 


ence system. These are scaled and plotted by a call to DRAW. 
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Pieumeswo—7 sana 5-8 demonstrate the plot 

options available to TARGET. 
e. Subroutine LAUNCH 

LAUNCH causes a plot display of the target 
ship's defensive missile envelope to be drawn. Two options 
are available. They are a three-dimensional cylinder or a 
bold black circle plotted on the surface as shown in Figures 
3-1 and 3-2, respectively. Selection of the desired 
envelope plot 1S a uSer input parameter. Inputs to LAUNCH 
are the missile range and the display type. 

The cylindrical display is defined using 
cylindrical coordinates such that the height of the 
cylinder cooresponds to twice the maximum aircraft input 
altitude. Cylinder construction consists of ten circles 
vertically separated by a distance of one fifth the maximum 
aircraft altitude. Each circle consists of 180 dots. 

The range ring display merely consists of 180 
coordinate points expressed in plane polar coordinates. The 
radius of the circle corresponds to the maximum missile 
launch range. 

Variables used in LAUNCH are: 

(1) RANGE - Maximum missile launch range. This 
is an input parameter from the calling program (SEAPLT). 

(2) IDISPL - Integer display option flag such 
that IDISPL=1 causes the single range ring display. Other- 


wise the cylinder display is drawn. 
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SER PLOT 


PERSPECTIVE VIEW COORD: 22708. -1389. 1600. 
€ACH RANGE AING «= sooo. METERS 
LAUNCH ENVELOPE QT 12900. METERS 





Paes o-7: IJARGET output. Ship heading=000 


SEA PLOT 


PERSPECTIVE VIEW COOAOs. 22708. -13689. 1600. 
CACH RANGE RING = S000. METEAS 
LAUNCH ENVELOPE AT 12000. METERS 


Figure 3-8: TARGET output. Ship heading=135 
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IV. SUMMARY AND CONCLUSIONS 


SEA PLOT is a three-dimensional computer graphics 
program written to display a war-at-sea encounter between a 
Single ship and a single aircraft, in which the ship is 
assumed to be the target and is located at the origin of an 
earth-fixed coordinate system. Roll, pitch, and yaw maneuvers 
of the attacking aircraft are shown along a three-dimensional 
flight path which depicts the attacker's coordinate positions 
in the earth-fixed reference system as well as the attitude 
of the aircraft body in three space. Up to 200 flight path 
coordinates can be plotted. To accomplish the projection of 
a three-dimensional scenario onto a two-dimensional plotting 
surface, a conic projection technique is used in which a 
Viewing plane and an eye position are automatically computed 
based upon the attacking aircraft's initial position. 

Displays available to the user include: (1) the earth- 
fixed coordinate axes, (2) 5000-meter circular range rings 
centered on the target, (3) ship defensive missile launch 
envelope displayed in solid cylindrical or plane circular 
Pomme 4) aircraft air track, (5) aircraft ground track, 
(6) aircraft maneuvers of roll, pitch, and yaw, (7) ship 
target body display located at the earth-fixed origin. To 
generate these displays, SEA PLOT was written using thirteen 


subroutines which fall into four functional categories: plot 
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management routines, graphics (vector drawing) routines, 
conic projection routines, and plotting (picture segment) 
routines. This organization was selected to allow for easy 
conversion from the present VERSATEC graphics software to 
other graphics languages such as the Tektronix developed 


rool 80. 
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I. INTRODUCTION 


This appendix discusses the Tektronix graphics software 
package PLOT 80 and its use on the Tektronix 4081 graphics 
computer. Because Reference 6 provides an adequate dis- 
cussion of the Tektronix 4081 equipment and system operation, 
computer light-off and operation will not be reiterated 
here. Instead, discussion will be focused on the graphics 
language PLOT 80 and its uses in the Tektronix 4081 program 
mode. 

PLOT 80 is a particularly powerful graphics tool con- 
Sisting of 374 different subroutines covering a wide range 
of graphics concepts including bit stream manipulation, 
vector drawing, test input and output, and picture display. 
These subroutines are discussed in detail in the PLOT 80 
user's manual, Reference 7. Because PLOT 80 represents a 
complex, sophisticated software system, Reference 7 should 
be reviewed by a prospective user prior to any attempt at 
programming of the Tektronix 4081. This appendix will 
discuss some of the more basic concepts employed by PLOT 80 
and provide the programmer some useful background informa- 
tion which is intended to complement the information in 
Reference 7. 

PLOT 80 will be discussed here in the context of the 


4081 program mode of operation as opposed to the host mode 


ad 





Meecwssed in Reference 6. In the program mode, all computa- 
tions are accomplished using the hardware organic to the 

4081 system itself; unlike the host mode in which a host 
computer, such as the IBM 3033 or PDP-11, performs all calcu- 
lations and uses the 4081 for graphics display only. The 
word size of the 4081 does inhibit many mathematical computa- 
Maomos such aS the transforms used in SEA PLOT, which result 
tmeunaerflows and hence limits to some extent the program 
mode applications of the 4081. However, this mode does 
provide an excellent tool for investigating graphics 
Structures such as the target and aircraft body representa- 
tions used in SEA PLOT. Section II of this appendix dis- 
cusses some of the important basic concepts of PLOT 80. 
ecemron JI] will describe the procedures for file creation, 


compiling, linking, and execution on the Tektronix 4081. 
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Elie RICE ex 


The 374 distributed graphics support subroutines (DGSS) 
in PLOT 80 are FORTRAN callable routines that permit the 
user to: 

(1) Perform sophisticated graphics in integer or 
floating point format. For example, routines that perform 
graphic transformations or create polygons, circles, and 
splines are included in the software package. 

(2) Perform graphic input (GIN) in windowspace 
Or viewspace coordinates. 

(3) Perform file manipulation. 

(4) Communicate with the host computer. 

(5) Perform assembly language operations including 
bit and character string manipulation, primitive I/O, 
and dynamic memory allocation, 

Before performing any DGSS graphics, a program must 
include the two following steps: 

CA BRREV (logical unit, device). 

CALL GONEW (logical unit). 

The first routine assigns a logical unit (to be defined 
Himeorlowmle paragraphs) to a graphic output device; 1l.e., 
the display controllery The» second routine creates a data 
Structure called the graphics control block (GCB) for the 


specified logical unit. By means of the logical unit, 


les, 





eiererore, a GCB is linked with a graphic output device as 


Shown in Figure A-l1. 


LOGICAL UNIT 







CORAL eC =OULPUT DEVICE 


Figure A-1l 


For example, if the programmer decides that the display 
ewmesorler 1S to be assigned to logical unit 0, he could 
create a GCB as follows: 

Cwnteoryey (0 CeMpC.') 
CALL GONEW (0) 

He is now ready to create the graphics control block for 
his particular graphics needs. 

The graphic control block contains a sequence of 
parameters that determine exactly where on the surface of 
the graphic output device the graphics information defined 
iMmanprogram is to be displayed. The GCB is essential to 
all DGSS graphics. Attempting to output graphics information 
to a device without first creating a GCB causes an error 


condition to be signaled which will terminate execution. 
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nl = Soe 


No DGSS graphics programmer can entirely ignore the 
concept of a GCB. The GONEW routine, however, allows a 
programmer to get directly to the problem of the graphics 
application without spending much time and space on 
creating a GCB and defining each and every one of its 
parameters. 

The GONEW routine automatically allocates the memory 
required for a GCB and initializes the contents to a set of 
default values. The programmer does not need to know where 
in memory the GCB is located or how much space it occupies. 
Table 14-1 of Reference 7 lists the GCB parameters and the 
default settings. 

Data transfers under the graphics operating system are 
always performed through a logical unit assigned to a 
device. The programmer communicates with the logical unit 
rather than the device itself. This process greatly simpli- 
fies the I/O operation. The programmer does not need to 
do "housekeeping" tasks, such as addressing the device and 
checking the status of the device to see if it is busy or 
ready to send and receive data. All housekeeping is per- 
formed by the proper device driver which is called into 
action when a device is assigned to a logical unit. There 
BaeomlGe1ogtcal Units, numbered 0 to 15. Initially, all units 
are assigned to the keyboard unit until a BFDEV call is made. 
Any logical unit may be assigned to any device. Graphic out- 


put devices are: 
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peedisplay control ler 
ee plotter 
(3) a file 


Another important concept available to the programmer 
is that of a picture segment. Using this concept, the over- 
all plotted output is divided into collections of vectors 
called picture segments. Each segment is numbered and can 
be specified independently of the other segments. The pro- 
grammer can therefore manipulate each segment independently 
of the others. (Recall that SEA PLOT was organized along 
these lines with the coordinate axes, surface grid, launch 
envelope, and flight path each being a separate picture 
segment capable of individual manipulation.) 

Picture segments are also sometimes called refresh 
objects, which refers to how the segment is displayed on the 
screen. Segments can be displayed in either refresh or 
storage. When a segment is drawn in storage, the segment 
is drawn once on the screen and remains visible until the 
screen is erased. A segment displayed in refresh, however, 
is not "stored'' on the screen. The Segment must be con- 
tinually drawn and redrawn on the screen to remain visible 
(a function automatically carried out and not a programmer 
mesponSibility) . 

Refresh graphics allow the modification and movement 
of segments without erasing the screen. The effect is 
yamiars co anamation in that a picture 1S drawn at one point 


on the screen, then redrawn again at another, then again at 
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another, and so on. The picture Seems to be moving across 
Pies creen. 
iijene are three steps in creating a picture segment: 
(1) Initialize the segment assigning it a number 
Mpeeco 255). 
(2) Specify vector attributes (coordinates, intensity, 
ete.) . 
(3) Terminate segment creation. 
These are accomplished by the OPEN and CLOSE subroutines, 
lemiliustrate the use of the GCB and picture Segments, 
a possible animation routine for SEA PLOT would be as follows: 
CNP Ee DEY (0,'DEC:*) 
CALL GONEW(0) 
C. DRAW COORDINATE AXES (PICTURE SEGMENT 1) 
CALL OPEN (1) 
CALL AXIS (1) 
gah CLOSE (1) 
C. PLOT AXES IN STORAGE 
GAmioee LX. (1) 
C. DRAW RANGE RINGS (PICTURE SEGMENT 2) 
CALL OPEN (2) 
GArL €GRID 
Gari CLOSE (2) 
C. PLOT RANGE RINGS IN STORAGE 


@Abh FIX (2) 
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C. DRAW AIRCRAFT (SEGMENT 3) IN REFRESH MOVING FROM 
C. MILESTONE TO MILESTONE 
DO 10 I=1,NPOINT 
READ X,Y,Z,ROLL, PITCH ,HEAD 
CALL OPEN (3) 
CALL ACRAFT(Z,Y,Z,ROLL,PITCH,HEAD) 
CALL CLOSE (3) 
C. PLOT SEGMENT 3 IN REFRESH 
CALL POST (3) 
C. TIME DELAY CORRESPONDING TO AIRCRAFT VELOCITY 
CALL TDELAY 
C. ERASE PREVIOUS BODY PLOT 
CALL UNPOST (3) 
C. ITERATE THROUGH REMAINING MILESTONES 
10 CONTINUE 
STOP 


END 
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Cie bee Renl LON AND EXEGUION 


This final section discusses Tektronix 4081 program mode 
procedures to create, compile, link-edit, and execute a 
program. Unlike the job control language options of the 
NPS IBM 3033 main frame computer, compiling, linking, and 
executing jobs on the 4081 are entirely distinct operations 


which must be individually performed by the operator. 


pee LE CREATION 

Creating a local file on the 4081 is accomplished by 
the command "EDIT filename. file type.'' This operation is 
essentially the same, using the same functions, as the EDIT 
system used on the main frame. When creating a multi- 
Structured file, one must first create a library file. 

This library file will then contain each of the individual 
subroutines comprising the overall program, Moreover, each 
subroutine must be individually compiled. 

Using SEA PLOT as an example, one would create a 
library file named SEAPLT containing all of the applicable 
Subroutines as follows: 

FORMAT SEAPLT.LIB(4) 

This creates a SEAPLT library with four data blocks 
available for use. 

To create the individual subroutine files in the 


SEAPLT library, the programmer must edit, or create, a 
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Seyeeate file in the library for each subroutine as indi- 


cated below: 
EDIT SEAPLT/MAIN. FOR 


INPUT 


ED ie orveb /CONPRO. FOR 
INPUT 


EDIT SEAPLT/TARGET.FOR 
INPUT 


Or so on until all subroutine files have been entered. 


B. COMPILING 


Each file structure in the library file must be 
individually compiled in the form: 
FORT library name/file name.OBJ, 
LST=library name/filename.FOR 
Compiler errors will be in the LST file, which can be 


read by typing the command 


Bend lol. LST 
Because this becomes somewhat tedious for a file 
with several structures, one can create a batch routine 


which will compile the various subroutines. The following 
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Socine will create a compile routine for a batch processor 
1 @ec 
EDIT COM. BAT 
LUBE 
FORT LIBRARY NAME/<1>.0OBJ, 
LST=LIBRARY NAME/<1>.FOR 
Exiy 
Using this routine, successive subroutines may be rather 
more expeditiously compiled by submitting the batch job 
BATCH COM. filename 
which will compile a single subroutine and must be repeated 


for each subroutine in the library. 


C. LINK-EDITING 
After a program is compiled, the program must be link- 
edited by typing the following commands: 
LINK 
OUT library name/RUN.OBJ 
LINK library name/main program name 
EDIT library name 
Ee to? DGL, SYS: RIL 
Ponta the Output of the link editor indicate that there 
are undefined terms in the program, this process must be 
repeated until all undefined terms have been resolved. This 
sometimes involves repeating the link-edit step four or 


five times, 
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Pee EXECUTION 
Executing the successfully compiled and linked program 
Mee retreshingly, a very straightforward procedure. One 
merely enters the command 
USR: library name/RUN 
If any execution errors are detected, they will be dis- 
played on the screen along with any graphics output which 


has been generated before the error was detected. 
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